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Abstract 

Using the classical estimation method of moments, we propose a new semiparametric estima- 
tion procedure for multi-parameter copula models. Consistency and asymptotic normality of 
the obtained estimators are established. By considering an Archimedean copula model, an 
extensive simulation study, comparing these estimators with the pseudo maximum likelihood, 
rho-inversion and tau-inversion ones, is carried out. We show that, with regards to the other 
methods, the moment based estimation is quick and simple to use with reasonable bias and 
root mean squared error. 
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1 Introduction 

Recently, considerable attention has been paid to the problem of inference about copulas. The 
monographs of Cherubini et al. (2004), Nelsen (2006) and Joe (1997) summarize to some extent the 
activities in this area. Roughly speaking, a copula function is a multivariate distribution function 
with uniform margins. It is used as a linking block between the joint distribution function (df) F 
of a vector of random variables X = (X\, X d ) and its marginal df's Fi, F d . This probabilistic 
interpretation of copulas is justified by the famous Sklar's theorem (Sklar, 1959) which states that, 
under some mild conditions, there exists a unique copula function C, such that 

F(x 1 ,...,x d ) = C(F 1 (x 1 ),..,F d (x d )). 
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In other words, the copula C is the joint df of the random vector U = (U\, ...,Ud) , with Uj = 
Fj (Xj) . That is, for u = (u±, Ud) , we have 

C(u) = F(Ff 1 (n 1 ),...,F- 1 M), 

where F~ l (s) := inf {x : Fj (x) > s} denotes the generalized inverse function (or the quantile 
function) of Fj. 

A parametric Archimedean copula model arises for X when the copula C belong to a class C := 
{Cq, 6 £ 0} , where O is an open subset of W for some integer r > 1. Statistical inference 
on the dependence parameter 6 is one of the main topics in multivariate statistical analysis. 
Several methods of copula parameter estimation have been developed, including the methods 
of concordance (Oakes, 1982, Genest, 1987), fully maximum likelihood (ML), pseudo maximum 
likelihood (PML) (Genest et al, 1995), inference function of margins (IFM) (Joe, 1997, 2005), 
and minimum distance (MD) (Tsukahara, 2005). The performance of the PML procedure vis-a-vis 
to the other methods has been discussed by several authors. For example, the simulation study 
carried out by Kim et al. (2007) has concluded that the PML method is conceptually almost the 
same as the IFM one. It overcomes its non robustness against misspecification of the marginal 
distributions. Moreover, by using the PML method, one would not lose any important statistical 
insights that would be gained by applying the IFM. An advantage of the PML over the IFM is that 
the former does not require modeling the marginal distributions explicitly. Therefore, the PML 
estimator is better than those of the ML and IFM in most practical situations. However, in time- 
consuming point of view the ML, PML, IFM and MD methods require intensive computations, 
notably when the copula dimension increases. Moreover, when using these methods the copula 
density has to be involved, therefore a serious inaccuracy at boundary points arises. Several 
numerical methods are proposed to solve this problem, but they are still inefficient when dealing 
with high dimensional copula models, more precisely for d > 2 (see, Yan, 2007, Section 5). 

The aim of this paper is to propose an alternative estimation method similar to the concordance 
one, avoiding technical problems caused by copula density and providing estimators with reason- 
able time-consuming, bias and root mean squared error (RMSE). The concordance method, also 
called the r-inversion and p-inversion, which are based, respectively, on Kendall's r and Spear- 
man's p rank correlation coefficients, used to estimate parametric copula models with at more two 
parameters. Indeed, the r-inversion and p-inversion methods use the functional representations 
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of r and p in terms of the underling copula C (Schmid et at, 2010), given by 



p = p(C) 



i){ 2d L c(u)<iu -| 



2 d - (d + 

More precisely, suppose that copula C is a parametric model, i.e. C = Co, then both r and p 
become functions in 6 as well, that is r = r (0) and p = p (6) . Let r and /5 be, respectively, 
empirical versions of r and p pertaining to the sample (Xi, X n ) from the random vector X and 
suppose that Cq is one-parameter copula model (i.e. r = 1). Then, the estimators of 6 obtained 
by r-inversion or p-inversion methods are defined by 6 := r _1 (t) or 9 := p _1 (p) , where r _1 and 
p^ 1 are the inverses, if they exist, of functions 9 — >■ r (9) and 9 — > p(6) respectively. In the case 
when r = 2, that is when 6 = (9\,9<2) , we have to use jointly the two inversion methods, called 
(t, ^-inversion, to have a system of two equations 

t(9 1 ,9 2 )=?, p(9 1 ,9 2 ) = p. (1) 

The consistency of such estimators is discussed in the Appendix section IA.3I In conclusion, when 
the dimension of parameter 6 equals r, we have to use r measures of association, for example 
Blomqvist's beta (3, Gini's gamma 7, ... (see, Nelsen, 2006, page, 207) which, in general, is 
not convenient on the choice of measures point of view. More precisely, suppose that we are 
dealing with a parameter 6 = (61,62) of a copula model Cq, then one has the right to ask the 
following question: What couple among all measures of association have to be chosen to get 
a better estimation for 01 On the other hand, it is worth mentioning that often there exist 
difficulties while using Spearman's rank correlation coefficient. One such difficulty is when using 
very large or very small samples. For example, in the case of very large samples, it is very time 
consuming to perform Spearman's coefficient since it requires ranking of the data of all variables. 
Then we have to look for an alternative more convenient class of measures providing estimators 
with nice properties. A solution to this problem may be given by applying the classical method 
of moments to random variable (rv) C (U) . Indeed, let us define the feth-moment (C) , called 
copula moment, of rv C (U) as the expectation of (C (U)) fc , that is 



M k (C) := E \(C (U)) k ] = f (C (u)) k dC (u) , k = 1, 2, ... (2) 
Notice that the case k = 1 corresponds to 

Mi(C) = E[C(U)] 



f2 d - 1 - 1) r + 1 
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In other words, M k (C) may be considered as a generalization of Kendall's rank correlation r. To 
our knowledge, the method of moments is only used in one-parameter copula models, also known 
by the r-inversion method (see for instance, Tsukahara, 2005). Note that, since < C (u) < 1, 
then Mk (C) are finite for every integer k. Now we are in position to present a new estimation 
method that we call copula moment (CM) estimation. Suppose that, for unknown parameter 
eO CW, we have C = Cg, then M k (C) = M k (0) , where 

M k (0) := [ (C e (u)) fc dC e (u) , k = 1, 2, ... (3) 

J[0,l] d 

From equations ([3]) , we may consider M : — > (Mi (0) , M r (0)) as a mapping from O CW 
to W\ that will be used as a means to estimate the parameter 0. More precisely, for a given 



sample (Xj, ...,X n ) of the random vector X, let us denote as the estimator of defined by 
( M k)i<k<v That is 

:= M (JMi, M r j , (4) 

where M k is the empirical version of M k (C) and M _1 is the inverse of the mapping M, provided 
that it exists. The rest of the paper is organized as follows. In Section [21 we present the main 
steps of the copula moment estimation procedure and establish the consistency and asymptotic 
normality of the proposed estimator. In Section [21 an application to multiparameter Archimedean 
copula models is given. In Section [H an extensive simulation study is carried out to evaluate and 
compare the CM based estimation with the PML and (r, /?)-inversion methods. Comments and 
conclusion are given in Section [5j The proofs are relegated to the Appendix. 

2 Copula Moments based estimation 

In this section we present a semiparametric estimation procedure for the copula models based on 
the CM's ([3]) . First suppose that the underlying copula C belongs to a parametric family Cq, with 
9 = (#i,--- ,9 r ), and satisfies the concordance ordering condition of copulas (see, Nelsen, 2006, 
page 135), that is: 

for every 6>i, 2 6 O : Q x + Q 2 ^ C 01 (> or <) C 02 . (5) 
It is clear that this condition implies the well-known identifiability condition of copulas: 

for every 9 U 2 G O : e 1 + Q 2 ^Ce 1 + Cg 2 . 
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Identifiability is a natural and even a necessary condition: if the parameter is not identifiable then 
consistent estimator cannot exist (see, e.g., van der Vaart, 1998, page 62). 

For a given sample (Xi, X n ) from random vector X = (X\, ...,X d ) , we define the corresponding 
joint empirical df by 

n 

F n (x) = rT 1 1 {X u < xi, ...,X di < x d } , 



i=l 



with x := (xi, ...,Xd) , and the marginal empirical df's pertaining to the sample (Xj\, .... Xj rt ) , 
from rv Xj, by 

n 

F jn (xj) = n" 1 ^2 1 {Xji < xj} ,j = l, d. (6) 



8=1 



According to Deheuvels (1979), the empirical copula function is defined by 

C n (u) := F n (F^ 1 (m) , F" 1 (u d )) , for u G [0, l] d , 

where F~^ (s) := inf {x : Fj n (x) > s} denotes the empirical quantile function pertaining to df 
Fjn- We are now in position to present, in three steps, the semiparametric CM-based estimation: 

• Step 1: For each j = l,...,d, compute Uji := Fj n (Xji) , then set 

Ui := (pu, Udi) , i = 1, n. 



Step 2: For each k = 1, ...,r, compute 

n 



i=l 



(7) 



as the natural estimators of CM's given in equation 



• Step 3; Solve the following system 



M 1 (6 1 ,...,6 r ) = M 1 
M 2 (e 1 ,...,9 r ) = M 2 

M r (9 1 ,...,e r ) = M r . 



(8) 



-CM 



The obtained solution := I 6\, 9 r ) is called the CM estimator for 0. 
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Consistency and asymptotic normality of Q C are stated in Theorem [T] below whose proof is 
relegated to the Appendix lA.il For convenience we set 

L k (u; 6) := (C e (u)) fc - M k (6) and L (u; 6) = (L t (u; 6) , L r (u; 6)) . (9) 

Let 6q be the true value of and assume that the following assumptions [H.l] — [H.3] hold. 

• [H.l] 6q 6 O C W is the unique zero of the mapping — > Jj Q L (u; 0) dCo (u) which is 
defined from to l r . 

• [H.2] L (•; 6) is differentiable with respect to 9 with the Jacobian matrix denoted by 



L (u; 0) is continuous both in u and 6, and the Euclidian norm 
dCg-integrable function h (u) . 



L (u; 6) 



is dominated by a 



• [-ff.3] The rxr matrix j4o := Jj jid L (u; Oq) dCg (u) is nonsingular. 

Theorem 1 Assume that the concordance ordering condition ([S]) and assumptions [H.l] — [H.3] 
hold. Then with probability tending to one as n — > oo, i/iere exists a solution to the system 
([8]) which converges to 6$. Moreover 

^(e OM -e ) %Ar(o,A^D (a^) t 



as n — > oo, 

w/iere D := var {L (£; O ) + V (£; O )} and V (£; O ) = (Vi (£; ) , V r (£; O )) with 

A r d(c eo ( ~" k 

pi i[o,i] d d Uj 

where £ := —,^) is a (0, l) d -uniform random vector with joint df Cq q 



V k (&0 ):=y2 [ a(C g' (u)) (l{^<u,}- Uj )dCg (u), k = l,...,r, 



Remark 1 The asymptotic variance A n 1 Dn (A n 1N ) T may be consistently estimated by the sample 



Ai--=J d L ( u ! J dCgow (u) and Di := L (u^; J+V(U;;0 

as is done, in Genest et al. (1995) and Tsukahara (2005) in the case of PML's estimator and 
Z- estimator respectively. For more details on the Z- estimation theory, one refers to van der Vaart 
(1998), page 41. 
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3 Application: Archimedean copula models 



As application to the CM estimation method, we consider the Archimedean copula family defined 
by C(u) = p~ l (^2j = i (p(uj)j , where tp : [0, 1] — > M is a twice differentiable function called the 
generator, satisfying: p (1) = 0, ip' (x) < 0, <p" (x) > for any x £ (0, 1) . The notation ip~ l stands 
for the inverse function of (p. Archimedean copulas are easy to construct and have nice properties. 
A variety of known copula families belong to this class, including the models of Gumbel, Clayton, 
Prank, ... (see, Table 4.1 in Nelsen, 2006, page 116). Let K c (s) := P (C (U) < s) , s G [0, 1] , be 
the df of rv C (U) , then equation (J2]) may be rewritten into: 

M k {C)= f s k dK c (s), k = 1,2, .... 
Jo 

Suppose now, for unknown 8 £ O, that ip = p e , it follows that C = Cq, Kc = K# and (C) = 
M k (6) , that is 

M k {6) = I s k dK e (s), k = 1,2,..., 
Jo 

Notice that, one of the nice properties of Archimedean copula is that the df Kc of C (U) may be 
represented in terms of the first and second derivatives of the generator. Indeed from Theorem 4.3.4 
in Nelsen (2006), for any s £ [0, 1] , Kg(s) = s — p e (s) jp' e (s) , it follows that the corresponding 
density is K' e (s) = ip'g (s) tp g (s) / (tp' e (s)) 2 . Therefore the kth CM, defined in (J2J> , may be rewritten 
into 



M k (6) 



1 kVe ( s ) i s ) 



ds, k = 1, 2, ... 



(10) 



'0 H(a)) 4 

In terms of Kg, the assumptions [H.l] — [H.3] and Theorem [1] may be rephrased, respectively, to 
[-ff.l'] — [-H"-3'] and Theorem [2] below. For convenience, we set 

£ (t; 0) = (£i (t; 6) , C r (t; 6)) with C k (t; 6) := t k - M k (0) . 

• [H.l'] 9q € O C W is the unique zero of the mapping — > j C (t; 6) dK.g (t) that is defined 
from O to W. 

• [-ff-2'] C (•; 0) is differentiable with respect to 6 with the Jacobian matrix denoted by 

~dM k (6) 



C(t;6) :-- 



89, 



C (t; 6) is continuous both in t and 9, and the Euclidian norm 
dKg-integrable function h (t) . 



C(t;0) 



is dominated by a 
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• [H.3'] The r x r matrix Aq := J C (t; 6q) dKg (t) is nonsingular. 

Theorem 2 Assume that concordance ordering condition © and assumptions [H.1'] — [H.3'] hold. 
Then with probability tending to one as n — > oo, i/iere exists a solution 6 to the system © 
which converges to Oq. Moreover 



as n oo, 



^ (e CM _ o ) 4 (o, ^-^o (Vf) , 

w/iere 

P := «ar j£ (£ 6> ) + ^ g (t) (1 {£ < i} - *) ^ (i) J , 

where £ is a (0, 1) -uniform rv and g(t) := _1 )i<t< *s r -dimensional vector. 

3.1 Illustrative example 

The Gumbel family is an Archimedean copula defined by 



C/3(u) = exp 



V 



^(-lnu^] |,/3>1, 



with generator tpg (t) = (—hity , (3 > 1. For the sake of flexibility in data modeling, it is better 
to use the multi-parameters copula models than the one-parameter ones. To have a copula with 
more than one parameter, we use, for instance, the transformed (or distorted) copula defined by 

Cr(u) = r- 1 (c(r(n 1 ),...,r(n d ))), 

where T : [0, 1] — > [0, 1] is a continuous, concave and strictly increasing function with T (0) = and 
T (1) = 1. As an example, suppose that V = T a , with T a (t) = exp (1 — t~ a ) , a > and consider 
the Gumbel copula Cp, then the transformed copula C aj p (u) = T^ 1 (Cp (T a (u\) , ...,T a (ud))) is 
given by 

d 



C a,p (u) := 



V 




;n) 



which is also a two-parameter Archimedean copula with generator <p a p (t) := (t~ a — 1)^ . Note 
that C a) p verifies the concordance ordering condition © (see, Nelsen, 2006, page, 145). By an 
elementary calculation we get the kth. CM: 

(fc+l) 2 ^ + (/t+l)Q'/3 



In particular the first two CM's are 

„ , 2/3 + a/3-1 , o . 3/3 + a/3 -2 

Ml C 01 ^ := 4/3 + 2a/3 ^ ^ := 9/3 + 3a/3 ' 

Let (Xi, X n ) be a sample of random vector X = (X\, ...,Xd) , then the CM estimator ( a, f3 
of (a, /3) is the unique solution of the system 

Mi (a, 0) = Mi 
M 2 (a,/3) = M 2 . 



That is 



^ 8Mi-9M 2 -l 3 l-4Mi+3M 2 

a = — — , P = t ~ — \ 7 ^rv- (12) 

l-4Mi+3M 2 (l-2Mi) (l-3M 2 



4 Simulation study 

First notice that all numerical computations are performed on a personal computer with a micro- 
processor speed of 2.4 GHz. To evaluate and compare the performance of CM's estimator with the 
PML and (r, /^-inversion estimators, a simulation study is carried out by considering the trans- 
formed bivariate Gumbel copula family C Qj/ g defined above. The evaluation of the performance is 
based on the bias and the RMSE defined as follows: 

N / , N 2 \ 1 / 2 



RMSE= l]T(^-tf) 2 . (13) 



Bias = . „ 

N ^— ' \N 

i=l \ i=l / 

where 0{ is an estimator (from the considered method) of from the ith samples for iV generated 
samples from the underlying copula. In both parts, we selected N = 1000. The procedure out- 
lined in Section [2] is repeated for different sample sizes n with n = 30, 50, 100, 200 to assess the 
improvement in the bias and RMSE of the estimators with increasing sample size. Furthermore, 
the simulation procedure is repeated for a large set of parameters of the true copula C Qj( g. For each 
sample, by using formulas (fT2"j) , we obtain the CM-estimator (Zxi,/3A of (a, 0) for i = 1, N, and 
the estimators a and j3 are given by a = i Y^hLi ®i an£ l P = ~k YliLi Pi- The choice of the true 
values of the parameter (a, /3) have to be meaningful, in the sense that each couple of parameters 
assigns a value of one of the dependence measure, that is weak, moderate and strong dependence. 
In other words, if we consider Kendall's r as a dependence measure, then we should select values 
for copula parameters that correspond to specified values of r by means of the equation 



r(a,/3) = 4/ C a p (m, u 2 ) dC a>/3 u 2 ) -1. (14) 
The selected values of the true parameters are summarized in Table 1: 
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T 


a 




0.01 


0.1 


1.059 


0.2 


0.2 


1.137 


0.5 


0.5 


1.600 


0.8 


0.9 


3.450 



Table 1: The true parameters of transformed Gumbel copula used for the simulation study. 







r = 


0.01 






r = 


0.5 






r = 


0.8 








a = 


0.1 


/3 = 


1.059 


a = 


0.5 


/3 = 


1.6 


a = 


0.9 


/3 = 


3.45 




n 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


CPU 


30 


-0.081 


0.330 


0.032 


0.180 


-0.051 


0.654 


0.039 


0.481 


-0.073 


0.907 


-0.372 


1.130 


22.013 sec 


50 


-0.046 


0.253 


0.022 


0.139 


-0.043 


0.487 


0.018 


0.367 


-0.032 


0.723 


0.261 


0.916 


49.563 sec 


100 


-0.026 


0.173 


0.009 


0.097 


-0.023 


0.350 


0.012 


0.262 


-0.027 


0.548 


-0.089 


0.733 


2.789 mins 


200 


-0.011 


0.122 


0.002 


0.064 


-0.009 


0.243 


0.006 


0.180 


0.003 


0.386 


-0.056 


0.506 


10.370 mins 


500 


-0.005 


0.075 


0.000 


0.041 


-0.007 


0.155 


0.003 


0.117 


-0.007 


0.241 


-0.026 


0.323 


1.035 hours 



Table 2: Bias and RMSE of CM estimator of two-parameter transformed Gumbel copula. 



a = 0.1 = 1.059 a = 0.2 = 1.137 a = 0.5 = 1.6 a = 0.9 = 3.45 

Bias RMSE Bias RMSE Bias RMSE Bias RMSE Bias RMSE Bias RMSE Bias RMSE Bias RMSE 

n = 30 

CM -0.082 0.313 0.039 0.190 -0.064 0.642 0.033 0.486 -0.068 0.567 0.074 0.466 -0.072 0.955 -0.359 1.128 

PML -0.067 0.068 -0.485 0.486 -0.117 0.128 -0.561 0.568 0.078 0.394 -0.426 0.594 -0.043 0.421 0.253 1.027 

p-T 1.236 2.895 -0.213 1.775 1.101 2.984 -0.913 1.641 -0.556 2.312 -0.975 1.142 -0.439 0.691 0.919 1.039 

n = 50 

CM -0.046 0.245 -0.021 0.141 -0.052 0.468 0.029 0.357 -0.037 0.506 0.015 0.364 -0.033 0.732 0.206 0.892 

PML -0.060 0.062 -0.478 0.482 0.102 0.112 -0.516 0.526 -0.072 0.240 -0.472 0.556 0.025 0.315 0.330 0.772 

p-T 1.115 2.033 -0.289 1.378 1.035 2.537 -0.354 1.341 -0.478 2.110 -0.952 1.021 -0.392 0.508 0.801 0.991 

n = 100 

CM -0.022 0.171 0.009 0.100 -0.019 0.342 0.016 0.258 0.029 0.367 -0.005 0.257 -0.025 0.551 0.155 0.704 

PML -0.058 0.059 -0.483 0.485 -0.109 0.113 -0.524 0.528 -0.071 0.158 -0.403 0.469 -0.023 0.167 -0.017 0.237 

P-t 0.973 1.220 -0.176 1.273 0.923 2.335 -0.340 1.457 -0.365 1.114 -0.852 1.001 -0.255 0.397 0.708 0.686 

n = 200 

CM -0.016 0.122 0.002 0.064 -0.011 0.244 0.001 0.383 -0.014 0.245 0.005 0.180 -0.001 0.396 -0.060 0.527 

PML -0.041 0.062 -0.503 0.505 -0.099 0.102 -0.514 0.516 -0.050 0.116 -0.333 0.415 -0.059 0.144 0.038 0.400 

p-T 0.874 1.025 -0.235 1.215 -0.890 2.414 -0.330 1.041 -0.321 0.997 -0.786 0.988 -0.239 0.331 0.580 0.625 



Table 3: Bias and RMSE of CM, PML and r-p estimators of two-parameter transformed Gumbel copula. 



5 Comments and conclusions 



From Table 2, we conclude that by considering three dependence cases: weak (r = 0.01) , moderate 
(r = 0.5) and strong (r = 0.8) , the performance, in terms of bias and RMSE, of the CM based 
estimation is well justified. In each case, for small and large samples, the bias and RMSE are 
sufficiently small. Moreover, in time-consuming point of view, we observe that for a sample size 
n = 30 and for N = 1000 replications, the central processing unit (CPU) time to process CM's 
method took 22.013 seconds, which is relatively small. For one replication N = 1, the CPU 
time (in seconds) for different sample sizes are summarized as follows: (n, CPU) = (30, 0.213) , 
(100, 0.312) , (200, 0.844) , (500, 3.922) . Table 3 shows that both the PML and the CM based 
estimation perform better than the (r, /^-inversion method. However, in weak dependence case 
r = 0.01, the CM method provides better results than the PML one, mainly when the sample 
size increases. On the other hand, it is worth mentioning that our method is quick with respect 
to the PML one. The main advantage of our method is that it provides estimators with explicit 
forms, as far as Archimedean copula models are concerned. This is not the case of the other 
methods which require numerical procedures leading to eventual problems in execution time and 
inaccuracy issues. In conclusion, the CM based estimation method performs well for the chosen 
model. Furthermore, its usefulness in the weak dependence case particularly makes it a good 
candidate for statistical tests of independence. 

Acknowledgement. The authors are indebted to an anonymous referee for valuable remarks 
and suggestions. 
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A Appendix 

A.l Proof of Theorem Q] 

By considering CM's estimator as a Z-estimator (van der Vaart, 1998, page 41), a straight ap- 
plication of Theorem 1 in Tsukahara (2005) leads to the consistency and asymptotic normality 
of the considered estimator. Indeed, the existence of a sequence of consistent roots to @ , 
may be verified by using similar arguments as the proof of Theorem 1 in Tsukahara (2005). More 
precisely, we have to check only the conditions in Theorem A. 10.2 in Bickel et al. (1993). Indeed, 
first recall ([9]) and set 

,. n 

$ (0) := / L (u; 0) dCg (u) , and $ n (0) := n^ 1 ^ L (u*; 

J fd i=l 

where Uj = (F\ n (Xu) , F^ n (X<&)) , with (Xji, ...,Xj n ) is a given random sample from the rv 
Xj. In view of assumption [H.2] the following derivatives exist 



1 



$(0)= / L(u;0)dC7 eo (u), $ n (0) = - V L U,; 
.hd n » V 



Next, we verify that 



sup 



$ n (0) - $ (0) 



: |0 — 0o| < e n > — )• 0, as n — > oo, 



(15) 



for any real sequence e n — > 0. Indeed, since L is continuous in 0, then 



sup 



and the fact that 



L Uf, — L XJf, 9, 



: |0 - O | < e n } = o P (1) , i = 1, n, 



$ n (0) - $ n (0 O ) 



1 



Tl 

n ^ 

i=l 



L[U,:;0] -L(U i; o 



implies 



sup 



— 0q | < e n r ~~ > 0' as n — ^ oo. 



$ n (0) - (0 O ) 

On the other hand, in view of the law of the large number, we have 

1 n • • 

- V L (U»; O ) 4 $ (0 O ) , as n -> oo, 



2 = 1 



(16) 



where Uj = {i<j (-Xji)}i<j<<2 • Moreover, in view of the continuity of function L in u and Glivenko- 
Cantelli theorem, that is 

sup \Fj n (xj) — Fj (xj)\ — > 0, j = 1, d, almost surely, as n — > oo, 
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we have 



n 



i=l 



It follows that (0 O ) - * (0o ) ^ 0, which together with (USD, implies ([15]). Conditions (MGO) 



and (MG3) in Theorem A. 10.2 in Bickel et al. (1993) are trivially satisfied by our assumptions 
[HI] — [H3] . In view of the general theorem for Z-estimators (see, van der Vaart and Wellner, 1996, 
Theorem 3.3.1), it remains to prove that \fn — (0o) converges in law to the appropriate 
limit. But this follows from Proposition 3 in Tsukahara (2005), which achieves the proof of 



A. 2 Proof Theorem H 

The proof of Theorem [2] is straightforward by using similar argument as the proof of Theorem [TJ 



A. 3 Consistency of (r, p)-inversion estimators 

In this section we give assumptions on copula models Ce, satisfying condition ([5]), that allow 
consistency of (r, /))-inversion estimators of = (9i,6 2 ) defined in ([TJ) . On the other terms, we 
propose some conditions of copula family Ce ensuring, for large sample sizes, both existence and 
uniqueness of system ([TJ . The idea is to express (r, /?)-inversion estimators in terms of Z-estimators 
(van der Vaart, 1998, page 41) and then we use similar assumptions allowing consistency of these 
estimators (see, van der Vaart and Wellner, 1996, Theorem 3.3.1). Indeed, recall that Kendall's 
tau and Spearman's rho corresponding to the couple of rv's (Xi,X 2 ) of dependence function Ce 
are defined, respectively, by 



Theorem [TJ 



□ 



therefore the details are omitted. 



□ 



p(0) = 12 u\UidCe {u\,ui) — 3. 

I J [o,i] 2 

It is easy to verify that (UJ) is equivalent to the following system 




n 



L T (F ln (X u ) , F 2n (X 2i ) ; 0) = 



(17) 



L p (F ln (X u ) , F 2n (X 2l ) ; 0) = 




L T (ui,u 2 ;0) := 4C (u 1 ,u 2 ) - 1 - r (0) , 
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and 



L p (ni, u 2 ; d) ■= 12uiu 2 — 3 — p(0) 



with Fj n denotes the empirical df pertaining to the sample (X/i, Xj n ) defined in §6§ . This 
implies, by representation (fTT|) , that (r, /^-inversion estimators of the true value 6q are, indeed, 
Z-estimators. Therefore, by using the Z-estimation theory, we conclude that consistency of such 
estimators may be established provided that the following two assumptions hold: 



[A.l] 6q, element of an open O C M 2 , is the unique zero of the mapping 



9 



low 



L (u 1 ,u 2 ;6) dCo Q (ui,u 2 ) , 



defined from O to R 2 , with L {u\, u 2 ; 0) := (L T (ui, u 2 ; 0) , L p {u\, n 2 ; 0)) . 

[A.2] L (•; 6) is differentiable with respect to with the Jacobian matrix denoted by 



L(ui,ti 2 ;0) := 



8L T (u\,U2\Q) dL T (u\,U2\0) 

89! ae 2 

dL p (ui,U2;9) dL p (ui,U2\9) 

ae 1 ae 2 



L(ui,u 2 ;0) is continuous both in (ui,ii 2 ) and 0, and the Euclidian norm 
dominated by a dCe-integrable function g (ni,u 2 ) . 



L(ui,u 2 ;0; 



is 
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